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Thirty years ago the Hubbard model was introduced by Gutzwiller, Hubbard and 
Kanamori with the main purpose of mimicking the ferromagnetism of transition 
metals. Soon after, Nagaoka and Thouless pointed out a basic mechanism for fer- 
romagnetism in strongly correlated electron systems by studying the motion of a 
single hole in a half-filled Hubbard model. This important work was hoped to shed 
light onto metallic ferromagnetism from the low doping regime. Unfortunately, this 
low doping route towards ferromagnetism has not been successful as far as rigor- 
ous results for finite doping concentrations are concerned. In the work presented 
here, we start from the opposite limit of low particle concentrations. In this limit 
we provide the first proof of a fully polarized metallic ground state for a Hubbard 
model. The proof proceeds by mapping Hubbard "zigzag" chains onto a continuum 
model with an additional degree of freedom and local first Hund's rule coupling. For 
this model the maximum total spin multiplet is shown to be the unique ground state 
for infinite Hubbard coupling. Our proof may open a low density route towards the 
understanding of the ferromagnetism of Hubbard models. 

PACS numbers: 71.27.+a, 75.10.Lp 

The Hubbard model was introduced about 30 years ag c&ii mainly for the pur- 
pose of mimicking the ferromagnetism of transition metals (see also the the paper 
by R. Strack and D. Vollhardt in this volume). In view of this it is quite aston- 
ishing that up to this point our understanding of the ferromagnetism of Hubbard 
models is very limited. All_thc rigorous results proving the presence of a ground 
state with net polarizatiortmQ are either concerned with half-filled bands (i.e. with 
ferromagnetism of insulators) or with somewhat artificial flat-band systems which 
have ferromagnetic ground states even for vanishing interaction strength. Very re- 
cently Tasaki succeeded in praying local stability of the ferromagnetic ground state 
for nearly-flat-band systems J3 There is not a single rigorous result which would 
provide,a, positive answer towards the basic question which was in the minds of the 
fathersoBB of the Hubbard model: 



• Is a strong intra-atomic Coulomb repulsion sufficient for pushing a partially 
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filled (i.e. paramagnetic) band into a fully polarized ferromagnetic ground 
state (Nagaoka state)? 

Nagaoka's theoremBi was, undoubtedly, an early important contribution to 
the subject. It tells us the precise conditions for the formation of ferromagnetic 
polarons around electrons inhales doped into a half-filled Hubbard model at infinite 
coupling strength-isee alsot2rliIrla). The puzzling results obtained for two and more 
electrons or holestJ made obvious, though, that the problem of the global alignment 
of ferromagnetic polarons, in case they existed, is much harder to solve. In fact, not 
a single rigorous result establishing a ferromagnetic ground state for finite doping 
concentration is known. The low doping route towards metallic ferromagnetism 
which was so impressively opened by Nagaoka's theorem has failed as far as rigorous 
results are concerned. 

There are, however, several pieces of evidence in favor or in disfavor of ferromag- 
netism for Hubbard models on various lattices. Let us first look at the hypercubic 
lattices of dimensions d = 1 to d = oo. For the linear chain (d = 1) it is exactly 
known that there is never a ferromagnetic ground state, except for the limit of in- 
finite coupling strength where due to complete spin degeneracy the fully polarized 
states are among the many ground statesllj In the limit of infinite dimension a 
variational argument using Gutzwiller wave functions alsOj-oroves that the ground 
state is never fully polarized, not even for infinite coupling.Efl Here, partially polar- 
ized ground states have not been excluded and probably exist, but have not been 
established so far because the numerical solution of the mean-field theory at d = oo 
(see the paper by A. Georges in this volume) is hard in the very strong coupling 
limit. Variational estimates of the regime of stability of a Nagaoka ground state 
suggest that the lower the dimension of the lattice the larger this regime is, if it 
exists. Therefore, this problem has been studied most extensively for the square 
lattice. The variational proofs of instability of the Nagaoka state for the square 
lattice leave a rather small regime of possible stability which can be characterized 
by a critical electron or hole concentration (relative to half-filling) of 5 cr = 0.29 
and by a critical r.ni.mlipg of U cr = 63t where t denotes the nearest neighbor hop- 
ping amplitude.li3li3'li3E3 Theuestimates of the ground state energy deduced from 
a high temperature expansions allow for a critical density 5 cr (again relative to 
half-filling) of at most 0.10 and are consistent with 5 cr =p-Q. This result appears not 
in agreement with finite square plaquette Lanczos dataE3 and with densitsi-matrix 
renormalization group (DMRG) extrapolations from multi-chain systemstj which 
both suggest a 5 cr of about 0.20. The evidence in the case of the square lattice seems 
still controversial and not conclusive. The even more favorable case of the double 
chain has no±_been studied as extensively as the square lattice. Lanczos datacj and 
DMRG dataEy agree on S cr = 0.22 for the double chain. It would be very helpful if 
a high temperature expansion for the double chain also became available. 

In the presence of so much uncertainty in the case of hypercubic lattices it 
should not be overlooked that there is positive evidence for Hubbard models being 
much more fond of ferromagnetism on non-bipartite lattices.lij : E3 To avoid confu- 
sion in view of the absence of particle-hole symmetry for non-bipartite lattices we 
assume the nearest neighbor hopping amplitude t to be positive (which implies the 
band minimum to be found at the center of the Brillouin zone). It is then demon- 
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strated with rather simple variational wave functions that for less than half-filling 
(n < 1) the Nagaoka state is never a ground state of thcJJubbard model. Fox 
more than half-filling (n > 1) both variational estimatesEj'Ej and Lanczos dataEHl 
suggest sizeable regimes with a Nagaoka ground state. This is perfectly consistent 
with the Nagaoka theorem which predicts ferromagnetic polarons for n > 1 only.cl 
The triangular lattice which is a prototype two-dimensional non-bipartite,lattice 
has been investigated with particular care. The best variational calculations^ esti- 
mate the critical density at infinite coupling U as n cr < 1.887, while Lanczos data 
appear to confirmEJ that there is a Nagaoka ground state for 1 X. jo. 1.76. For 
the three-dimensional fee and hep lattices similar investigational^' E3'E3 suggest a 
Nagaoka ground state for the whole regime 1 < n < 2 at sufficiently large coupling. 
The same type of behavior follows from Lanczos data for a quasi-one-dimensional 
lattice, the chain of triangles.t3 We thus have good evidence for the rather striking 
phenomenon of fully polarized ground states in the limit of small hole densities 
(relative to complete filling, i.e. to n = 2). It is this phenomenon which we are 
going to use in this paper as a starting point for opening a new route towards the 
understanding of the ferromagnetism of Hubbard models. Since one can always 
map the limit of small hole densities (n — > 2) onto the limit of small particle den- 
sities (n — > 0), employing a particle-hole transformation, we will call this route the 
"low density route" in contrast to what we called the "low doping route" in the 
discussion above. 

In what follows we will consider the t — t' Hubbard chain 

n = -T,{ tc ^+^ +t ' c l^ + h - c ) + u T, n ^- w 

ia i 

For the special choice t' = t this model describes the chain of triangles mentioned 
above. We do prefer, in fact, to visualize this system as a zigzag chain of atoms 
rather than as a linear chain because we will have to satisfy the condition t' > t/A 
to obtain a ferromagnetic ground state. This condition can probably only be met if 
the second neighbor distance is shorter than twice the first neighbor distance. With 
this in mind we also avoid any conceptual conflict with the "no-go-theorem" for 
ferromagnetism in strictly one-dimensional systems.Efl 

As we want to focus on the low density limit of the zigzag chain we imagine that 
a particle-hole transformation and a gauge transformation have been performed 
such that from now on we will be interested in the case t > and t' < 0. Looking 
at the dispersion 

e(fc) = -2t cos k - 2t' cos 2k (2) 

of the t — t' band we observe a qualitative change at the bottom of the band as 
t' decreases: the band minimum at k — splits into two equivalent minima for 
t' < — 1/4. We will see that in the low density limit the existence of these two band 
minima is responsible for the formation of a Nagaoka ground state in the strong 
coupling limit. 

The two minima form around momenta ±ko where cosfco = —t/At'. At low 
densities only states in the neighborhood of these two minima matter. It is im- 
portant for the following argument that due to the ultra-violet stability of one- 
dimensional systems which implies intrinsic high energy cut-offs this is also true for 
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the interacting system. In the presence of two band minima the system therefore 
develops an internal degree of freedom (left valley - right valley) at low densi- 
ties. Introducing the notation c\. ra := c£_ rfco a (t = ±1) and fitting cosine bands 
e T (k) = I - 2tcos(k - rk ) with e = 2t' + t 2 /4t' and t = -4t'[l - {t/At') 2 } into the 
two minima of the original dispersion (Q) we express the original Hamiltonian TL in 
terms of new fields distinguished by the symbol ~. The result of this simple change 
of notation will be discussed using the new Hamiltonian 



where Si T and fii T denote the spin and particle density operators of state r = ± at 
site i. 

The Hamiltonian TL describes two Hubbard chains with nearest neighbor hop- 
ping t which are coupled to each other by the on-site exchange interaction J. It has 
twice as many degrees of freedom as the original Hamiltonian ([!]). Nevertheless, TL 
and TL are closely related in the limit of small particle densities: for J = 2U the 
two models have one and the same continuum limit. The ferromagnetic exchange 
interaction of strength 2U results from the inter-valley scattering contained in the 
Hubbard term of the original model. 

The presence of a ferromagnetic exchange interaction J > in TL suggests to 
expect ferromagnetic ground states for this model. It is, in fact, possible to arrive 
at the following rigorous statement: 

• At U — oo, for any J > and for any particle number < N < 2L all the 
ground states of TL have maximum total spin S = N/2. 

Here L is the number of sites per chain such that N = 2L corresponds to half- 
filling of the double chain system. The proof of the above statement is quite simple. 
Let us first consider the_case J = where we have two uncoupled Hubbard chains. 
From the exact solutionEa we know that at U — oo the ground states are given by a 
Slater determinant of spinless particles for each chain combined with a completely 
arbitrary spin wave function. Since the kinetic energy is a convex function of the 
particle number the particles distribute equally among the two chains (N + = N- 
for even N and N + = 7V_ ± 1 for odd N). Including now the exchange term we 
see that the energy of the ground states is increased unless particles are in a triplet 
state whenever they meet at the same site. The only way to make sure that particles 
never meet in a singlet state is to have maximum total spin. This simple argument 
is easily made rigorous by a formal proof that out of the subspace of completely 
spin degenerate ground states of the uncoupled chains only those with a totally 
symmetric spin wave function are uneffected by the positive semi-definite exchange 
term. 

The proof that TL has no other than fully polarized ground states is an inter- 
esting result in its own right. Interpreting the internal degree of freedom as an 
orbital degeneracy we can rephrase this result by saying that a Hubbard chain with 





(3) 
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orbital degeneracy and first Hund's rule coupling has Nagaoka ground states for 
U = oo. We are reminded of the fact that orbital degeneracy combined with ferro- 
magnetic Hund's rule coupling is a genuine feature of real ferromagnets in the first 
place. Probably, omitting orbitaLdegeneracy in modeling metallic ferromagnets by 
the single band Hubbard modeltlBu was an oversimplification after all. In cases 
where the Hubbard interaction alone is not sufficient to drive a system ferromag- 
netic Hund's rule coupling might serve to provide the final kick. 

Despite of this remark we now turn back to the original single band model. Due 
to the equivalence of models discussed above we conclude that the continuum limit 
of H has only fully polarized ground states at t' < — 1/4 and U = oo. This result 
represents the first proof of metallic ferromagnetism in a Hubbard model. Applying 
a scale transformation to the continuum model we find the following formula for 
the ground state energy per site at given densities of particles n and magnetization 
m: 

E /L = (e - 2t)n + n 3 t$(U/nt, m/n). (4) 

We don't know the scaling function fi) since the model is not integrable. But 
from the above proof we infer 

$(oo,/x) >*(oo,l) (p<l). (5) 

The above result opens a low density route towards the understanding of the 
ferromagnetism of Hubbard models on non-bipartite lattices. We have learned 
that in the presence of several valleys at the bottom (or top, after a particle- 
hole transformation) of the band the Hubbard interaction implies a ferromagnetic 
exchange interaction between particles in the various valleys. In one dimension this 
turned out to be sufficient for making the strong coupling system ferromagnetic. It 
is presently not obvious how this result can be generalized. Two valleys are also 
found ip-fhe band structure of the triangular lattice. The discussion presented above 
tells us£3 that in this case two valleys are not sufficient for ferromagnetism at low 
densities. Nevertheless, multi-valley situations seem to be a generic prerequisite for 
Hubbard models having ferromagnetic ground states in a broad range of particle 
densities. One may wonder whether the low density route towards ferromagnetism 
will turn out to be more successful than the low doping route. 
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